Coupled-channel model for charmonium levels 
and an option for X(3872) 
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The effects of charmed meson loops on the spectrum of charmonium are considered, 
with special attention paid to the levels above open-charm threshold. It is found that 
the coupling to charmed mesons generates a structure at the DD* threshold in the 
1 ++ partial wave. The implications for the nature of the X(3872) state are discussed. 

PACS numbers: 12.39.-x,13.20.Gd,13.25.Gv,14.40.Gx 

I. INTRODUCTION 

The charmonium spectroscopy has again become a very interesting field. On one hand, 
the 2 1 Sq i]' c state was discovered by Belle O and confirmed by BABAR [2] and CLEO (^J, 
and the l x Pi h c was observed in Fermilab |4| and by CLEO [5]. The masses of these long 
missing states are in perfect agreement with the predictions of quark model. At the same 
time, a new state A (3872) was found by Belle Q and CDF Q. 

This discovery has attracted much attention. As the state is just at the DD* threshold, 



it was immediately suggested 



that it might be a DD* molecule bound by pion exchange 



deuson"), considered long ago in [l| and, much earlier, in jl(3] and jllj]. This req uires 1 



12. 



++ 
The 



quantum numbers, and this assignment seems to be favoured by the data 
discovery channel is 7r + 7r~ J ftp with dipion most probably originating from the p. Together 
with the observation of X in the ujJ/ip channel 15j this opens fascinating possibilities for 
strong isospin violation, which is also along the lines of the deuson model. 



Other options for X(3872) are under discussion in the literature, see e.g. 
The most obvious possibility of X being a cc state seems to be ruled out by its mass: the 
state is too high to be a ID charmonium, and too low to be a 2P one [17| . This assumes 
that we do know the spectrum of higher charmonia, namely, the fine splittings and the role 
of coupling to Z)-meson pairs. It is the latter issue which is addressed in the present paper. 
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The mechanism of open-flavour strong decay is not well-understood. The simplest model 
for light-quark pair creation is the so-called 3 P model, suggested many years ago [^(j]. It 
assumes that the pair is created with vacuum ( 3 Pn) quantum numbers uniformly in space. 
The application of this model has a long history 21, 22, 13]. Systematic studies 0,13 of 
the decays of light and strange quarkonia show that with a 3 Pq -type amplitude calculated 
widths agree with data to within 25 — 40%. Recently the charmonia decays [26 1 and decays 
of D- and .D^-mesons J27I] were considered in the framework of the 3 Pq model. 

There exist also microscopic models of strong decays, which relate the pair-creation in- 
teraction to the interaction responsible for the formation of the spectrum, by constructing 
the current-current interaction due to confining force and one-gluon exchange. Among these 
is the Cornell model [3] which assumes that confinement has Lorentz vector nature. The 
model assumes that the confining interaction is the scalar one, while one-gluon-exchange 
is, of course, Lorentz vector. Possible mechanisms of strong decays were studied in the frame- 
work of Field Correlator Method (FCM) j^] , and an effective 3 Po operator for open-flavour 
decay has emerged from this study, with the strength computed in terms of FCM parameters 
(string tension and gluonic correlation length). 

Most of the above-mentioned papers are devoted to computing the widths, and only a 
few consider the effects of virtual hadronic loops on the spectra. The Cornell model [28[ has 
presented a detailed analysis of charmonia with coupling to D-mesons taken into account. 
The recent update j^l] of the Cornell model has presented splittings caused by coupling 
to mesonic channels for ID and 2P cc levels, confirming the previous result: A(3872) is 
well above the range of ID levels and well below the range of IP ones. The paper [32] has 
reported first results for hadronic shifts of lower charmonia due to mixing with D-meson 

argc. 
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pairs, calculated within the Pq model. The shifts appear to be alarmingh 
Meanwhile, phenomenological coupled- channel models like I33I I34I I39 
experience on the possibilities to generate nontrivial effects due to the coupling to hadronic 
channels. As a recent example one should mention the analyses [37] of new D s j states 
with masses considerably lower than quark model predictions, and coupling to mesonic 
channels being responsible for these anomalously low masses. It is interesting to note that 
;he coupled-channel calculations performed in the framework of chiral Lagrangian approach 
I has arrived at the same conclusions. 

In this paper the coupled-channel model for charmonia levels is presented, based on 



3 



the nonrelativistic quark model for cc spectrum and 3 Po _ type model for pair-creation. In 
Section II the dynamics of coupled channels is briefly outlined. Section III introduces the 
quark model. Sections IV and V contain the results which are discussed in Section VI. The 
paper ends with a short summary. 



II. DYNAMICS OF COUPLED CHANNELS 

The details of coupled-channel model can be found e.g. in [3], [2^. Here I review the 
essentials. 

In what follows the simplest version of coupled-channel model is employed. Namely, it is 
assumed that the hadronic state is represented as 

i*)=L.^i^,,,y (i) 



.EiXi|Afi(i)M 2 (i)> 

where the index a labels bare confined states \ip a ) with the probability amplitude c a , and 
Xi is the wave function in the i-th two-meson channel \M 1 (i)M 2 (i)). The wave function 
obeys the equation 

f H c V \ 

H\$) =M|tf), H= . A , (2) 

\ V H MlM2 J 

where H c defines the discrete spectrum of bare states, with H c \tp a ) = M a \ip a ). The part 
Hm\M2 includes only the free-meson Hamiltonian, so that the direct meson-meson interaction 
(e.g., due to t- or w-channel exchange forces) is neglected. The term V is responsible for 
dressing of the bare states. 

Consider one bare state |-0o) (the generalization to multi-level case is straightforward). 
The interaction part is given by the transition form factor /i(p), 

(MV\M iX M i2 ) = / t ( Pi ), (3) 

where Pj is the relative momentum in i-th mesonic channel. Then (J2J) leads to the system 
of coupled equations for c (M) and X«,Af(Pi) : 

co(M)M„ + yj f /i(p) Xi ,M(p)d 3 p = Mc (M), 
ma + m i2 + — J Xi,MKPi) + c o( M )/i(Pi) = MxiMVi)- 
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Here u,- = m ' im ' 2 is the reduced mass in the system of mesons with the masses m,-i and 
m i2 , and M is the mass of the bare state. In what follows the formalism will be applied 
to the system of charmed mesons, so the nonrelativistic kinematics is employed in (J3J). The 
fully relativistic version of coupled- channel model is presented in j^j]. 

With the help of (jlj) one easily calculates the t- matrix in the mesonic system: 

where E{ = M — inn — ma- The quantity g(M) is often called the hadronic shift of the bare 
state \ipo), as the masses of physical states are defined, in accordance with eq.(0), from the 
equation 

M - M + g(M) = 0. (6) 

Let the eq.(JH} have the solution Mb with Mb smaller than the lowest mesonic threshold, 
so there is a bound state with the wave function 

= ( ■ A (* fl |* fl > = l, cose = (^ B ). (7) 



\ f^ + Q) 2 / \ dM 



sin9J2i XwiPi)). 

Here ^XiMPi) is normalized to unity, and cos 9 defines the admixture of the bare state 
\ijjo) in the physical state \^b)- The explicit expression for this admixture reads 

-i / . . \ -i 

(8) 

-U r I - / \ i.J.\l I 

£i = ^1+171^2— M B , €i > 0. As far as I know, this Z-factor was first introduced by S.Weinberg 
in j^J many years ago as the field renormalization factor which defines the probability to 
find the physical deuteron \d) in a bare elementary-particle state \do), Z = |(<io|<i)| 2 . 

Even more detailed information is contained in the continuum counterpart of the factor 
Z, the spectral density w(M) of the bare state, given by 



w 



(M) = J2wi(M), Wi(M) = 47r/i^|c(M)| 2 e(M - m a - m i2 ), (9) 



where c(M) is the probability amplitude to find the bare state in the continuum wave 
function \^/)m- With c(M) found from the system of equations (jlj), one can calculate w(M): 

W<yM ^ = 2m (m-M + ^(M) ~ M-M + g(M)) ' (10) 



As shown in |4l| . the normalization condition for the distribution w(M) follows from the 
completeness relation for the total wave function (0) projected onto bare state channel, and 
reads: 

POD 

w(M)dM =1- Z, (11) 



'moi+mo2 

if the system possesses a bound state, and 



c OO 



/ w{M)dM = 1, (12) 

Jm 01 +m 2 

if there is no bound state (moi + is the lowest threshold). In the case of bound state 
present, all the information on the factor Z is encoded, due to eq. ([lip, in the w(M) too. 
On the other hand, the analysis in terms of w(M) can be performed in the case of resonance 
as well, as exemplified in |42j. 

In the latter case the t-matrix poles are situated in the complex plane. While the posi- 
tions of the poles are the fundamental quantities, another quantities are useful for practical 
purposes. Namely, one defines the visible resonance mass Mr from the equation 

M R - M + Re g(M R ) = 0, (13) 

and calculates the visible width as 

d Re g(M) 



T = Im g(M R ), 9fc = 1 



dM 



(14) 



M=M R/ 

Clearly this brings the t-matrix into Breit-Wigner form, i.e. in the form in which experimen- 
tal data are usually delivered. In what follows the factor 3ft will be called the renormalization 
factor. 

There are some limitations of course, as not the every peak has the Breit-Wigner shape. 
In the case of overlapping resonances the formulae ()13|) and (|14|) do not work. The special 
case of near-threshold S'-wave resonance is not described by Breit-Wigner or Flatte formula, 
and the scattering length parametrization is more appropriate j^ . 

The quantities Z and w(M) are the ones of immediate relevance. Indeed, there is no 
hope that, say, the elastic DD scattering will be measured some time. Our knowledge on 
mesonic resonances comes from external reactions, like e + e~ annihilation, 77 collisions, B- 
meson decays etc. etc. Assuming that such reactions proceed via intermediate qq states, 
one obtains that the cross-section is proportional to w(M): 

<t(— » mesons) oc r 0r w(M), (15) 
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where r 0r is the width of the bare state corresponding to the external reaction. Such formulae 
were used in j^gj to describe the e + e~ annihilation into charmed mesons. In the limit of 
narrow resonance eq. p5j) is reduced to the standard Breit-Wigner formula 

mesons) oc ^ (M _ ^ + ip , ^ 

where T = 2lmg(Mo) is the (small) width of the resonance. Similarly, for the bound state 
case the width T r for a given reaction is renormalized as 

T r = Z r°. (17) 

III. THE QUARK MODEL 

This section specifies the form factors fi(p)- The pair-creation model employed is the 
3 Po one, that is the pair-creation Hamiltonian is the nonrelativistic reduction of 

H q = g q J d 3 x^ q , (18) 

or a given flavour q, but two important points make it different from the model used in 

_ _ _ 

The approach |2J, |25|, |26| assumes that the pair creation is flavour- independent, which 



yields for the constant g q the form 

g q = 7 2m g , (19) 

where 7 is the effective strength of pair-creation. The factor 2m q implies enhancement of 
strange quarks creation comparing to light quarks one. There are no fundamental reasons 
to have such enhancement. Moreover, such factor is absent in microscopical models of pair 
creation, like 



281 ] and [29]. So, throughout the present study, I use the effective strength 7 
for the creation of light (u- and d-) flavours, while for strange quarks the effective strength 
7s = ^7 is used, where m q and m s are the constituent masses of light and strange quarks 
correspondingly. 
The authors of 



241 [25! . 13] argue that the assumption of flavour-independence gives a 
reasonably accurate description of known decays. One should have in mind, however, that 
the calculations (2^ 13 Q] are performed with the so-called SHO wavefunctions, i.e. with 
the wave functions of harmonic oscillator, and with the same oscillator parameter (3 for all 
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states. This assumption looks implausible, as the behaviour of form factors is defined by 
scales of wavefunctions, which, in turn, are defined by quark model. 

In what follows the standard nonrelativistic potential model is introduced, with the 
Hamiltonian 

2 . 

H = — + V{r) + C, V(r) = ar--^, (20) 
m c 3 r 

m c is the mass of charmed quark. This Hamiltonian should be supplied by Fermi-Breit- 
type relativistic corrections, including spin-spin, spin-orbit and tensor force, which cause 
splittings in the 2S+1 Lj multiplets. In the first approximation these splitting should be 
calculated as perturbations, using the eigenfunctions of the Hamiltonian (|20|) . The same 
interaction V(r) should be used in spectra and wavefunction calculations of D (D s ) mesons. 

In the first approximation the pair-creation amplitude is to be calculated with the eigen- 
functions of the zero-order Hamiltonian. Use of the SHO wavefunctions simplify these 
calculations drastically. So the procedure adopted is to find the SHO wavefunctions (of the 
form exp {—\(3 2 r 2 ) multiplied by appropriate polynomials) for each orbital momentum L 
and radial quantum number n, with the effective value of oscillator parameter (3 for each L 
and n, which reproduce the r.m.s. of the states. 

I use the following set of potential model parameters: 



0.55, a = 0.175 Gel/ 2 , m r = 1.7 GeV C = -0.271 GeV. 



m q = 0.33 GeV, m s = 0.5 GeV. (21) 
The spin- dependent force is taken in the form 

V sd = V HF + ^ L ■ S - L ■ S + i^L T, (22) 
m ^ r s lm A c r m A c r 6 

where L-S and T are spin-orbit and tensor operators correspondingly, and Vhf is the contact 
hyperfine interaction, 



^(r) = 6{r) S q ■ S 9 , (23) 



32na x ~ 
9m 2 . 

where, fC.owing the hnes of W , Gau.sian-smea.ing of the hyperfine interaction is intro 
duced, 

*(r) = (4=) e " KV ' ( 24 ) 

with k = 1.45 GeV. 
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TABLE I: Masses and effective values (5 (in units GeV). 



TlL 


B(nL) 


kJ Lj CI' Lj 


Mass 


is 


0.676 


l 3 5i 


3.264 








3.135 


25 


0.485 


2 3 5'i 


3.905 






^ '-'O 


3.850 


IP 


0.514 


l 3 Po 


3.773 






l 3 P 
1 J 1 


3.718 






l 3 Pn 


3.631 






1 X P 


3.732 


2P 


0.435 


2 3 P? 


4.230 






2 3 P 

J. 


4.181 






2 3 Pn 

J. (J 


4.108 






2 x Pi 


4.192 


ID 


0.461 


1 3 L» 3 


4.051 






1 3 £> 2 


4.043 






1 3 P>1 


4.026 






l 1 ^ 


4.043 



The masses and effective values of oscillator parameters (3 for the model ()21)1 are listed in 
the Table I. The effective values of oscillator parameter for .D-mesons are (3d = 0.385 GeV, 
and (3 Da = 0.448 GeV. 

One should not take the numbers given in last column too seriously, especially for higher 
states, as the fine splittings are not well-known, and the expression (}2*2*|) is surely too naive. 
Moreover, various much more sophisticated approaches, which reproduce the splittings in 



44|. 



IP multiplet, give different predictions for higher multiplets, as discussed in detail in 

The D-meson masses taken are M D = 1.867 GeV, M D * = 2.008 GeV, M Ds = 1.969 GeV, 
Md* = 2.112 GeV, so that the mass difference between neutral and charged .D-mesons is 
not taken into account. The pair-creation strength for light quarks 7 = 0.322 is used. The 
3 P amplitudes are listed in the Appendix A. 
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IV. LOWER CHARMONIA AND ^-LEVELS 

The single-level version of the coupled-channel model is used in what follows, with the 
exception of 2 3 Si — l 3 Di levels. 

All the physical charmonium masses below threshold are known. The hadronic shifts and 
bare masses were calculated from the equation 

M = M ph y S + 5, 5 = g(M phys ) = 9i{M phys ), (25) 

i 

where the sum is over mesonic channels DD, DD*, D*D*, D S D S , D S D* and D*D* (an 
obvious shorthand notation is used here and in what follows: DD* = DD* + DD*, and 

d s d* s = n.jy ■ njn. 

Besides, as the position of the 1 ^(3770) state is well-established, the bare mass of 
l 3 Di state was reconstructed by means of eq IT3*l and the visible width was calculated as 

©• 

The parameters of the underlying quark model (J2T|) and the pair-creation strength 7 = 
0.322 were chosen to reproduce, with reasonable accuracy, the model masses of IS, IP and 
2S states, and the width of ^(3770). 

The results for bound states are given in Table II together with corresponding values of 
Z-factors. The shifts are much smaller than in 2(|, but still substantial. 

Let me now discuss the 1 3 D\ level, lying above DD threshold. The bare mass is calculated 
to be 4.018 GeV. The calculated width of the ^(3770) is 25.5 MeV, which compares well 
with the PDG value of 23.6 ± 2.7 MeV j^J. Note that it is the visible width, while naive 
calculations would give 

T = 2 Im g(M R ) = 34.3MeV. (26) 

So the effect of coupling to mesonic channels on the width of ?/>(3770) is not small, 5i = 0.743. 

The mass of ^(3770) is less than 100 MeV higher than the mass of ^'(3686), with DD 
threshold opening in between, so one could in principle expect that these states are mixed 
due to coupling to mesonic channels. This mixing is not large in the given model, as 
the mixing disappears if the mass difference between various D-mesons is neglected (see 
the corresponding spin-orbit recoupling coefficients listed in Appendix A). The relevant 
formulae for two-level mixing scheme are given in the Appendix B, and here I quote the 
results. With physical masses of ^'(3686) and ^(3770) the bare masses are reconstructed 
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TABLE II: Hadronic shifts (in units MeV) of charmonium states below DD threshold due to 
individual channels, total shifts 5, bare masses and Z-factors. The results for the 2 3 5i state do 
not take into account mixing with the l 3 Z?i state 



n 2S+1 L i 


DD 


DD* 


D*D* 


D V D V 

-*— ' S S 


D.D* 


D*D* 

.s ^ s 


§ 


Mn 

■ LV - L \J 


z 


l 3 S 1 


11 


42 


69 


5 


19 


31 


177 


3274 


0.899 


1% 





59 


55 





26 


25 


165 


3145 


0.913 


1 3 P 2 


25 


64 


82 


8 


22 


27 


228 


3784 


0.804 


1 3 P 1 





70 


91 





22 


32 


215 


3726 


0.820 


l 3 ^o 


29 





118 


8 





43 


198 


3613 


0.841 







87 


76 





29 


27 


219 


3744 


0.817 


2 3 5i 


21 


60 


87 


45 


16 


27 


216 


3902 


0.743 


2^0 





83 


71 





24 


22 


200 


3838 


0.802 



as M(2 3 S 1 ) = 3.899 GeV, and M(l 3 £>i) = 4.016 GeV, so the masses are shifted only by a 
few MeV due to the mixing. The width of the ^(3770) becomes only 18.4 MeV, but it is 
the visible width, and the deviation from the value 25.5 MeV obtained without mixing is 
mainly due to the illegitimate attempt to fit the system of two overlapping states with a 
single-Breit-Wigner lineshape. The states indeed do overlap, as shown at Fig.l, where the 
spectral densities of bare 2 3 S' 1 and 1 3 -Di bare states are plotted. The lineshape of l 3 _Dx is 
distorted due to the mixing, and the lineshape of the 2 3 S'i is drastically changed, displaying, 
instead of small smooth background, a peak at the mass of about 3760 MeV. The 10 MeV 
difference between the peak position and the mass of ■0(3770) is again due to the prescription 
of visible width: single-Breit-Wigner approximation is not appropriate both for 3 Si and 3 D\ 
lineshapes. As to the lower state ^'(3686), the admixture Z\ of the 2 3 S'i bare state is 0.742, 
while the admixture Z 2 of the bare l 3 Di state is 0.00063. 

The calculated mass of bare 1 3 D 1 state is 4.016 - 4.018 GeV. So the value of 4.026 GeV 
given in the Table I looks quite acceptable. No other ID-states are known, so in what follows 
I take, with reservations mentioned in the previous section, the values of bare masses from 
the Table I. With the bare l 1 /^ state having mass of 4.043 GeV, the physical state is a 
bound state, 

M(1 1 D 2 ) = 3.800 GeV, Z(1 1 D 2 ) = 0.712. (27) 
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a) b) 

FIG. 1: a) Spectral density of 2 3 5i bare state with (solid line) and without (dashed line) mixing, 
b) The same for the bare l 3 Z?i state. 

Similarly, 

M(rlD 2 ) = 3.806 GeV, Z(1 3 D 2 ) = 0.689 (28) 

for the mass 4.043 GeV of the bare 1 3 P>2 state. The l 3 D% state is allowed to decay into DD, 
but the width is extremely small, as the DD system is in the P-wave: 

M(1 3 D 3 ) = 3.812 GeV, T(1 3 D 3 ) « 0.7 MeV, K(1 3 P> 3 ) = 0.717. (29) 

V. 2P-LEVELS 

Coupled-channel effects do not cause dramatic changes for the charmonia states discussed 
in the previous section. The situation with 2P-levels promises more, as 2P-charmonia are 
expected to populate the mass range of 3.90—4.00 GeV, where more charmed meson channels 
start to open, and some of these channels are the S-wave ones. 

The importance of S- wave channels follows from the eq.(jS}- The form factor f(p) of 
the S'-wave mesonic channel behaves as some constant at small p, so the derivative of the 
hadronic shift g(M) with respect to the mass is large for the masses close to the S'-wave 
threshold. As the result, hadronic shift due to the coupling to S-wave channel displays 
rather vivid cusp-like near-threshold behaviour. 

The physical masses and widths of 2P states were calculated with bare masses given by 
Table I, and the results are listed in Table III. Two different values are given for the mass 
and width of the 2 3 Pq state, for two different choices of the bare mass, see below. Looking at 
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TABLE III: The masses and widths (in units of MeV), and renormalization factors K of 2P char- 
monium states. The results for the 2 3 Po state are given for two different values of bare mass as 
described in the text 



State 


Bare mass 


Mass 


Mode 


Width 


K 


2 1 P 1 


4200 


3980 


DD* 


50 


0.615 


2 3 P 2 


4230 


3990 


DD 


22 


0.603 








D S D S 


1 










DD* 


45 










total 


68 




2 3 Pi 


4180 


3990 


DD* 


27 


0.543 


2 3 P 


4108 


3918 


DD 


7 


0.622 


2 3 P 


4140 


3937 


DD 


8 


0.393 



the numbers one would say that nothing dramatic has happened due to S-wave thresholds. 
Indeed, all the shifts are about 200 MeV, and the renormalization factors are about 0.5 — 0.6. 
It is the behaviour of spectral density which reveals the role of S'-wave thresholds. 

For the 1 l P\ case the S'-wave thresholds are DD* with the spin-orbit recoupling coefficient 
C = -1/V2 and multiplicity 4, the D S D* with C = -1/V2 and multiplicity 2, the D*D* 
with C — 1 and multiplicity 2, and D*D* S with C — 1 and multiplicity I. If the resonance is 
in the mass range 3.90 — 4.00 GeV (and it appears to be so), then the relevant thresholds are 
DD* and D*D*. The spectral density of the bare 2 1 Pi state is shown at Fig. 2. It displays 
the relatively steep rise near DD* threshold, and a beautiful well-pronounced cusp due to 
the opening of D*D* channel. 

In the 2 1 Pi case the S-wave strength is shared equally between DD* and D*D* channels, 
while for the 2 3 P 2 case all S-wave strength is concentrated in the D*D*, channel. As the 
result, the cusp due to the opening of D*D* channel is more spectacular in the 2 3 P 2 case, 
as shown at Fig. 3. 

The case of 2 3 Pi is even more interesting. Here, similarly to the 2 3 P 2 case, all the S-wave 
strength is concentrated in two channels, DD*, D S D*, and the multiplicity of the former 
is 4. So the strongest S-wave threshold is the DD* one, well below the resonance. The 
behaviour of the 2 3 Pi bare state spectral density is shown at Fig. 4, and is very peculiar: 
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w 




3.9 3.925 3.95 3.975 4 4.025 4.05 

FIG. 2: Spectral density of the bare 2 1 Pi state. 




FIG. 3: Spectral density of the bare 2 P2 state. 

together with a clean and relatively narrow resonance, there is a near-threshold peak, rising 
at the flat background. The DD* scattering length appears to be negative and large, 

%B« = -8 fm, (30) 

signalling the presence of virtual state very close to the DD* threshold, with the energy 
e = 0.32 MeV. So the coupling to mesonic channels has generated not only the resonance, 
but, in addition, a virtual state very close to physical region. 

The near-threshold peak should fade with the increase of bare state mass, and strengthen 
otherwise. There are uncertainties in the fine splitting estimates, so the mass of the bare 
2 3 Pi state could easily be about 10 MeV larger or smaller. The dependence of spectral 
density behaviour on the bare mass is shown at Fig. 5. The peak becomes less pronounced 
for the bare mass of 4.190 GeV, but the scattering length remains rather large, a —5.2 fm, 
which corresponds to the energy of virtual state of about 0.76 MeV (compare this with the 
scattering length in the l + ~ channel, \a\ « 1 fm). 10 MeV decrease of the bare state mass 
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FIG. 4: Spectral density of the bare 2 3 Pi state. 




FIG. 5: Near-threshold 2 3 Pi spectral density for the mass of bare state Mq = 4.180 GeV (solid 
line), M = 4.190 GeV (long-dashed line), M = 4.170 GeV (dashed line). 

leads to incredibly large scattering length, a ~ —17.8 fm, and virtual state with the energy 
0.065 MeV. Further decrease of bare state mass leads to moving the state to the physical 
sheet, i.e. to appearance of the bound state. This happens at the bare mass of about 4.160 
GeV, which seems, in the present model, to be beyond acceptable range for fine splitting. 
Similarly, the bound state appears if the pair-creation strength is increased by several per 
cent. 

The 2 3 P level is a disaster, as it always happens with scalars. The bare mass is con- 
siderably lower that the c.o.g., and the uncertainty in the fine splitting estimate is large. 
The S'-wave D*D* and D*D* channels are too high. The S-wave DD channel is too low. 
The only relevant S'-wave channel is D S D S (C = ^/|, multiplicity 1), and corresponding 
threshold is at 3.938 GeV, i.e. around the region where the resonance is expected. So, 
depending on the position of the bare state, variety of spectral density behaviour can be 
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FIG. 6: Spectral density of the bare 2 3 P 

state for the mass of bare state Mq = 4.110 GeV (solid 

line), M = 4.140 GeV (dashed line). 

achieved, as shown at Fig. 6, where spectral density is plotted for M (2 3 P ) = 4.108 GeV and 
Mo(2 3 Po) — 4.140 GeV. Note that two curves of Fig. 6 correspond to two completely different 
situations. The curve for Mq(2 3 Pq) = 4.108 GeV displays the normal resonance behaviour 
with a tiny cusp due to the opening of the D S D S channel. The curve for M (2 3 P ) — 4.140 
GeV is not resonance-like at all. The formal exercise of calculating the 9ft-factor and visible 
width does not mak e much sense, and, a 8 suggested in Q, sueh elation e„ r ve should be 
analysed in terms of scattering length approximation, and not in terms of Breit-Wigner or 
Flatte distributions. 



VI. DISCUSSION 

The quark model (J2T| is not the result of the sophisticated fit, it is rather a representative 
example. A serious fit should include proper treatment of fine and hyperfine splittings, as 
well as the mixing of bare 2 3 Si and l 3 Pi states. The calculations should be performed with 
more realistic wavef unctions, and not with the SHO ones. Relativistic corrections should 
be taken into account in calculations of bare spectra, and more realistic model should be 
used for P-meson wavefunctions. More QCD-motivated model should be employed for 
the pair-creation Hamiltonian, and loop integrals g(M) are to be treated relativistically. 
Nevertheless, there are gross features which are model- independent. 

In various pair-creation models, the shifts within the each orbital multiplet are approxi- 
mately the same, and differ only due to the mass difference of bare states in the multiplet 
and different masses of charmed mesons. This is model-independent, as the charmed quark 
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is heavy, and the wavefunctions within each nL multiplet do not differ much. The same is 
true for the wavefunctions of D and D* (and D s and D*) mesons due to heavy quark spin 
symmetry. 

In addition, the shifts for all states are more or less the same. There are no a priori 
reasons for this, but the numbers given in Table II are stable up to overall multiplier 7 as 
soon as the scales (3 behave as expected from quark model, and the value of 7 is constrained 
by experimental value of the ?/>(3770) width. 



Indeed, the present results appear to be very similar to the ones given in [28J, |3l|. One 
might suggest that, from phenomeno logical point of view, the effect of coupling to mesonic 
channels can be approximated by adding a negative constant to the potential. But it is 
not the whole story, as the coupling to charmed mesons generates the admixture of D- 
mesons in the wavefunctions of bound states, which affects such quantities as e + e~ widths 
in accordance with ec|17|) and the rates of radiative transitions, as discussed in detail in 



The e + e~ widths of J/ip and ^'(3686) are more or less accurately described by Van Royen- 
Weisskopf formula with QCD correction, so the renormalization (fTTj) is not harmless, even 
if it is as mild as 10 % reduction required for J/ip by the results of Table II. For ?/>'(3686) 
Z-factor is 0.743, so renormalization is larger. 

In the nonrelativistic quark model the e + e~ width of the 3 Di state is zero. The mixing 
communicates the e + e~ width of the 2 3 S\ state to the ^(3770) region, with the e + e~ — > DD 
crossection given by eqlTHI where the spectral density of 2 3 Si state ws(M) is substituted. 
It is reasonable to estimate the e + e~ width of the ■0(3770) using the peak value of Ws(M), 
which yields less than 1/3 of the measured value 0.26 ± 0.04 keV |45[. Similar result was 



obtained in 



The bare l 3 Di state has the small e + e width of its own as relativistic 



correction, and the bare 2 3 Si and l 3 Di states are to be mixed by tensor force. The scale 
of the admixture required to reproduce the relatively large e + e~ width of the ^(3770) is 
not small, as shown in 46]. The coupling to charmed mesons is able to explain only about 
1/3 of the observed e + e~ width of ^(3770), so direct mixing between bare 2 3 Si and l 3 Di 
bare levels is still needed, which would reduce the e + e~ width of the ifi' further. While the 
problem of leptonic widths is an open problem for coupled-channel model, it is clear that 
the values of Z-factors considerably smaller than given in Table II would destroy fragile 
agreement with the data on e + e~ widths achieved by quark model practitioners. 
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The 2 3 S' 1 -1 3 P 1 mixing due to coupling to meson channels is small, and this is also model- 
independent, as it vanishes when the masses and wavefunctions of P-mesons are taken to be 
the same. Note, however, that both ?/>(3770) and ip' states contain considerable admixture 
of four-quark component of their own, in the form of various P-meson pairs. Thus, if the 
mass difference between charged and neutral P-mesons is taken into account, some isospin 
violation could be generated. It is argued in that a small admixture of / = 1 component 
is needed to explain some discrepancies in the observed properties of "0(3770) and ip r . The 
question of whether such admixture can be generated by coupled-channel effects certainly 
deserves attention. 

The shifts of P-levels are more or less the same, the relevant thresholds are P-wave ones, 
and the 2 3 Si-l 3 Pi mixing does not affect the shift of the bare 3 Pi level drastically. So 
the combined effect of fine splitting and splitting due to coupled-channel effects on other 
P-levels is not large. In particular, it means that, with the mass of ^(3770) as experimental 
input, the physical 3 D 2 or 3 P 3 state cannot be placed as high as 3.872 GeV and, therefore, 
cannot be identified as X(3872), unless something drastic happens with fine splittings in the 
lP-multiplet. 

Due to the presence of S'-wave thresholds, the situation with 2P-levels is more interesting. 
The coupling to D*D* channels generates pronounced cusps in the spectral densities of bare 
2 l Pi and 2 3 P 2 levels, and the coupling to PP* channel generates the strong threshold effect 
in the 1 ++ wave. Within the given model, it is a virtual state with the energy less than 
1 MeV. The mechanism of generating such a state is quite peculiar: it is one and a same 
bare 2 3 Pi state, which gives rise both to the 1 ++ resonance with the mass of about 3990 
Mev, and a virtual state at the DD* threshold, i.e. where the X(3872) is observed. Due to 
the presence of strong S'-wave threshold, the hadronic shift appears to be large enough to 
destroy the one-to-one correspondence between bare and physical states, as widely discussed 
in connection with light scalar mesons 

To what extent this prediction is robust? Changes of the underlying quark model param- 
eters or of the value of pair-creation strength can shift this extra state either to the physical 
sheet or away from the physical region. In the latter case, however, the DD* scattering 
length remains large. One should have in mind that if such dynamical generation of extra 
state at DD* threshold is possible in the charmonia, the 1 ++ channel is the most appropriate 
place for this phenomenon. The latter statement is model independent. 
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First, note that the scalar charmonium does not decay into DD* at all, and the tensor 
one decays into DD* in the D-w&ve. As to 1 ++ and l +_ levels, they both have the desired 
5-wave decay mode. Apply now the heavy quark spin selection rule (3], which suggests 
that the spin of a heavy quark pair is conserved in the decay. The S'-wave decay mode comes 
from the four-quark state ccqq with all relative angular momenta equal to zero. Then the 
total angular momentum J = 1 can result only from the quark spins. The combination of 
S C c = 1 and Sgg = 1 is C-even, while the combination of S ce = and S gq - = 1 is C-odd. It is 
a simple algebra exercise to show that, symbolically, 

(cc) s= i ® (qq) s =i - ^( DD * + DD*), (31) 

and 

(cc) s=0 ® (qq)s=i - -\{D*D - DD*) + ^=D*D*, (32) 

and, independently of the pair-creation model, all the S'-wave strength of the 3 Pi decay is 
concentrated in the DD* channel, while in the l P\ decay it is shared equally between DD* 
and D*D* . Thus the threshold attraction in the 1 ++ DD* channel is always much stronger 
than in the 1 H one. It would be interesting to see if the pair-creation model of [2^ is able 
to generate lar ge D P* scattering length. 

The model [49] contains the detailed analysis of the X(3872) as a state bound both 
by pion exchange and quark exchange in the form of transitions DD* — > J/ip p,J/ip 
with the latter contributions being important for the binding. In fact, this model has 
predicted the J/ip u decay mode of the X(3872), and, after observation of the decay 
X(3872) — > 7t + 7t~7t°J/iI), it has become almost official model of X(3872). This is challenged 



by preliminary data from Belle 
7T + 7T~J/ip one, while the model 



14j on large D°D°7r° rate, more than ten times larger than 
49] claims the opposite. 
One could question validity of the naive quark-exchange model. Besides, the pion ex- 
change is definitely attractive in the 1 ++ channel, but there are uncertainties in the actual 
calculations; the details of binding depend on the cutoff scale A, as recognized in 4^. The 
attraction found in the coupled-channel model is large, and could help binding without large 
quark-exchange kernels, and, correspondingly, without large 7i + 7i^J/ip rate. 

From practical point of view, the wavefunctions of both models are not very distinguish- 
able. Indeed, the near-threshold virtual state of the coupled-channel model owes its existence 
to the bare 1 ++ state, but the near-threshold admixture of the bare state in the wavefunc- 
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tion is extremely small, as seen from Fig. 4 (recall that the spectral density is normalized to 
unity). So the decays like DDn and DD~f would proceed via D* decays, as described in [50]. 
As to short-distance decays and exclusive production, the rates of these are governed by 
large scattering length. This phenomenon was called low-energy universality in j^jj. Con- 
sider, for example, the near-threshold production of DD* pairs in the reaction B — > DD*K. 
As explained in Section II, the DD* invariant mass distribution is proportional to w(M) in 
the coupled channel model, with the lineshape plotted at Fig. 5. Now compare these curves 
with the ones 

amplitude, and observe that the low-energy universality indeed takes place. 

The lineshape for DD* production depends only on the modulus of scattering length, 
and the cases of bound state and virtual state are not distinguishable. If some inelastic 
channel like ^ is present, the, as shown in S, the lineshape for this channel depend, 
on whether there is a bound state or virtual state (cusp). The area under the cusp is 
much smaller than the area under the resonance, so, in principle, the data could distinguish 
between A(3872) as a bound state or a virtual state. 

To conclude the discussion of coupled-channel virtual state as A (3872) I note that, with 
7 MeV difference between D°D*° and D + D*~ thresholds, the coupled-channel model is able 
to generate considerable isospin violation. The coupled- channel calculations which resolve 
the mass difference between charged and neutral D-mesons are in progress now. 

Last year yet another new state was reported [3] as an enhancement in uiJ/ip mode, 
with the mass of 3941 MeV and the width of about 90 MeV (V(3940). As there are two 
1 ++ physical states per one bare 2 3 Pi, it is tempting to identify the 1 ++ resonance with the 
V(3940) state, following the suggestion of Nevertheless, with the given set of quark 
model parameters the resonance is 40 — 50 MeV higher than 3940 MeV, and is narrow. 

One might suggest the ++ assignment for V(3940), as the coupled-channel model places 
the scalar just at the right place. However, as the D S D S channel is opening at 3938 MeV, 
the width of the state cannot be as large as observed 90 MeV independently of what channel 
saturates it, unless, for some strange reasons, the state couples weakly to D S D S . It is similar 
to a (980)// (980) case, where, in spite of large coupling to light pseudoscalars (nn or %ri), 
the visible width is small due to strong affinity to KK threshold. 
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VII. CONCLUSIONS 

The spectrum of charmonium states below 4 GeV is calculated taking into account cou- 
pling to the pairs of lowest D and D s mesons. The analysis is performed within the 3 Pq 
model for light-quark pair-creation. It appears that quite moderate modification of quark 
model parameters is needed to describe charmonia below DD threshold, while coupling to 
open charm does not cause drastic effects on ID-levels. This is in contrast to 2P-levels, 
where opening of strong S'-wave channels leads to pronounced threshold effects. 

In particular, coupling of the bare 2 3 Pi state to DD* channel generates, together with the 
1 ++ resonance with the mass of 3990 MeV, a near-threshold virtual state with the energy 
of about 0.3 MeV, which corresponds to the extremely large DD* scattering length a »s —8 
fm, with the possibility to identify this state with X(3872). The admixture of the bare cc 
state in the near-threshold wavefunction is very small, so it is essentially the DD* state. 
The 1 ++ channel appears to be the only one where such state can be formed. 
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Appendix A 



The 3 Pq form factors are defined following the lines of [24j,|2j|, adapted for charmonium 
transitions. Let the cc meson A decay to qc meson B and cq meson C. Then the spin-space 
part of the amplitude in the cm. frame of the initial meson A is given by 

/(p) = / rf 3 ^(k - p) < s q s q \6(k)\0 > <^(k - r q p)<j>* c (k - r q p), (A.l) 

4>a is the wavefunction of the initial meson in the momentum space, cbs and (be are the 
wavef unctions of the final mesons B and C, p = p R , r„ = — ^ — , m c is the mass of charmed 

if f tji q m q +m c ' L 

quark, m q is the mass of light quark. O(k) is the 3 P operator: 

6(k) = -2 7 a ■ k. (A.2) 

The form factors for the transition A BC are given below in the IS basis, where / 
is the orbital momentum in the final meson system, and S = Jb + Jc is the total spin of 
the mesons B and C. In the narrow width approximation these form factors fig define the 
partial widths Ta^bc as 

T a ^bc = 2 Im g B c{M A ) = "^PbcHbc ^2 \hs\ 2 - ( A - 3 ) 

is 

The form factors calculated with SHO wavefunctions take the form 

fe ^ exp (-fSf)^ (A ' 4) 

where (3a and (3b = Pc are the oscillator parameters of initial meson A and final mesons B 
and C, 7 is the pair-creation strength, and r q = m q /(m c + m q ), m c is the mass of charmed 
quark, m q is the mass of light quark. The polynomial Vis is a channel-dependent one: 

Vis = fi C ls , (A.5) 
where Cis are spin-orbit recoupling coefficients for specific mesonic channels, and fi are: 
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-P 
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2 5/2 ^3 
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5/3 



Here 



/3| + 2r^l 
/3I + 2/31 



2 ' 



and 



P 2 



3/3l/3| 



(A.14) 



(A.15) 



+ 2/31 

For the transitions to strange mesons, one should replace r q by r s = m s /(m c + m s ), and 
insert the multiplier — , as explained in the main tex;. 

The decays of light mesons were considered in ^24\. which corresponds to r q — 1/2. In 
the case of /3a = Br = 3, and r q = 1/2 the expressions for the amplitudes listed above are 



equal to those of 



24j up to the factor of 1/2. In general, there are two graphs with different 
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TABLE IV: The sums \ Cis\ 2 for given initial and final states 



State 
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D*D* 






I = L - 1 


I = L + 1 


I = L — 1 


I = L+ 1 


I = L - 1 


/ = L + 1 






i 




9 




7 


Do 




n 
u 




q 
o 




D 


1 1 


n 


n 


1/2 


5/3 


1 


10/3 

-L\J j (J 


3 Po 
r 2 


n 


1 


n 


3/2 


2 


8/3 


r \ 


n 


n 


1 


5/fi 


n 


e; 

O 


3 Po 


3/2 











1/2 


20/3 


l D 2 








1/4 


7/5 


1/2 


14/5 







1 





4/3 
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3 D 2 








3/8 


14/15 


1/4 


56/15 


3 A 


5/12 





5/24 





1/6 


28/5 



topologies which contribute to the 3 Po amplitude, and the sum of both is quoted in [24], 
while in actual calculations each graph contributes with the individual flavour factor. In 
the case of charmonia transitions only one graph contributes, so that the amplitude for the 
transition into given charge channel is equal to (jA.4|) with the flavour factor of unity. 

The mass difference between neutral and charged mesons is not taken into account, so 
the sum over charge states is equivalent to introducing the multiplicity factor 2 for DD and 
D*D* channels, and 4 for DD* channel. The multiplicity factor for D S D S and D*D* is 1, 
and it is 2 for D S D* channel. 

Spin-orbit recoupling coefficients are tabulated in the Appendix A of [24J]. In the single- 
level coupled-channel calculations the squares of spin-orbit recoupling coefficients are needed, 
given in Table IV. Note that the multiplicity factors for DD* channels are twice as large as 
of DD and D*D* ones. 

The case of 2 3 Si-l 3 -Di mixing requires explicit expressions for spin-orbit recoupling co- 
efficients, as relative signs are important in the two-level mixing scheme. These are given 
below: 

Cio( 3 Si — s* 1 Sq + 1 So) = 1 
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C 10 ( 3 5 1 ^ 3 Si + 3 Si) = Jl 
Cn(( 3 S 1 ^ 3 ft + 3 Si) = -J— (A.16) 
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Cio( 3 -Di Sq + 1 So) - 

Cn{ 3 Di ^ 3 S 1 + 1 S ) = 

C W { 3 D X ^ 3 5! + 3 S{}-. 
C 12 ( 3 D 1 ^ 3 5j + 3 S 1 

C 32 ( 3 D 1 ^ 3 ^ + 3 Si) = -\/^ (A.17) 
Appendix B 

The formulae necessary to describe the 2 3 Si~l 3 Di mixing are collected here. 

Two sets of form factors, fs and fp, are introduced, which describe the transitions 
between mesons and 2 3 Si, l 3 Di levels. The system of coupled channel equation similar to 
(JU) leads to the DD t-matrix: 

t(p,p',M) =X)/^s(P)vW/^fl(P'). (B-l) 

where sum is over S and D states, and 

M -M D + g DD {M) 
Tss{M) = A(M) ' 

M - M s + gss(M) 



tdd(M) 



A(M) 



r DS (M) = r SD (M) = - 9 ~f^, (B.2) 

A(M) = [M - M 5 + <? SS (M)][M - M D + to2J (M)] - g 2 SD , (B.3) 
Ms and Md are the masses of bare states, and 

^(M) = 5>^(M), ft^(M) = / ^^ A, (B.4) 



27 



the index i labels mesonic channels. 

The visible physical masses are defined, similarly to ([13)1 . from the equation 

Re A{M phys ) = 0. (B.5) 

There are two solutions of this equation, M b below DD threshold, corresponding to the 
bound state, and M a above threshold corresponding to the resonance. The visible width is 
defined for the latter as 

r = 25? Im A(M a ), 
3T 1 = [1 + d ss (M a )][M a -M D + Reg DD (M a )] 



+ [1 + d DD (M a )][M a -M s + Reg S s(M a )} - 2Reg SD (M A )d SD (M A ) J 

d Re g^M) 
^ > dM 

The probabilities to find bare states in the wavefunction of bound state are 



Z s = [M b -M D + g DD (M b )]Z-\ Z D = [M b - M s + g ss (M b )]Z-\ (B.7) 

and Z~ l is obtained from 3? _1 given in (|B.6j) with the replacement M a — > M b . The spectral 
densities of bare states are: 

1 ( M-M D + g* DD M-Md±9dd\ m R , 

WsW = 2^\ A*(M) A(M) > (R8) 



and 



(M\ 1 f M-M s + g* ss M-M s + g S s \ m Q , 

^ (M) = 2^l A*(M) A(M) J' (R9) 



